Exact Analytical Solution of the Constrained Statistical Multifragmentation Model 
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A novel powerful mathematical method is presented, which allows us to find an analytical solution 
of a simplified version of the statistical multifragmentation model with the restriction that the largest 
fragment size cannot exceed the finite volume of the system. A complete analysis of the isobaric 
partition singularities is done for finite system volumes. The finite size effects for large fragments 
and the role of metastable (unstable) states are discussed. 
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I. INTRODUCTION 

Exactly solvable models with phase transitions play a 
special role in the statistical physics - they are the bench- 
marks of our understanding of critical phenomena that 
occur in more complicated substances. They are our the- 
oretical laboratories, where we can study the most funda- 
mental problems of critical phenomena which cannot be 
studied elsewhere. A great deal of progress was recently 
achieved in our understanding of the multifragmentation 
phenomenon [l|,ll>l3 when an exact analytical solution of 
a simplified version of the statistical multifragmentation 
model (SMM) H ||] was found in Refs. 6, 7]. This ex- 
act solution allowed us to elucidate the role of the Fisher 
exponent t on the properties of (tri)critical point and to 
show explicitly ^] that in the SMM the relations between 
T and other critical indices differ from the correspond- 
ing relations of a well known Fisher droplet model [2|. 
Note that these questions in principle cannot be clari- 
fied either within the widely used mean-filed approach or 
numerically. 

Despite this success, the application of the exact ana- 
lytical solution [a, Q| to the description of experimental 
data is very limited because this solution corresponds to 
an infinite system volume. To extend the formalism for 
finite volumes it is necessary to account for the finite 
size and geometrical shape of the largest fragments, if 
they are comparable with the system volume. It is clear 
that these corrections may be important for not too di- 
lute systems. Therefore, to have a more realistic model 
it is necessary to abandon the arbitrary size of largest 
fragment and consider the constrained SMM (CSMM) in 
which the largest fragment size is explicitly related to the 
volume V of the system. This will allow us to solve the 
CSMM analytically at finite volumes, to consider how the 
first order phase transition develops from the singulari- 
ties of the isobaric partition lO"! in thermodynamic limit, 
and to study the finite size effects for large fragments. 



II. LAPLACE-FOURIER TRANSFORMATION 

The system states in the SMM are specified by the 
multiplicity sets {uk} (rik — 0, 1, 2, ...) of fc-nucleon frag- 



ments. The partition function of a single fragment with 
k nucleons is 0: V(j)k{T) = V {■mTk/2TTf^'^ zu , where 
k — 1,2, ...,v4 {A is the total number of nucleons in the 
system), V and T are, respectively, the volume and the 
temperature of the system, m is the nucleon mass. The 
first two factors on the right hand side (r.h.s.) of the sin- 
gle fragment partition originate from the non-relativistic 
thermal motion and the last factor, z^, represents the 
intrinsic partition function of the k- nucleon fragment. 
Therefore, the function 4'k{T) is a phase space density 
of the k- nucleon fragment. For fc = 1 (nucleon) we take 
zi — A [A internal spin-isospin states) and for fragments 
with fc > 1 we use the expression motivated by the liquid 
drop model (see details in Ref. m): z^ = exp(— /fe/T), 
with fragment free energy 
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'W{T) k + a{T) fc2/3 + {t + 3/2)Tlnfc , (1) 



with W{T) = Wo + rVeo- Here Wo = 16 MeV is the 
bulk binding energy per nucleon. T'^ je^ is the contri- 
bution of the excited states taken in the Fermi-gas ap- 
proximation (eo — 16 MeV). o[^^ is the temperature 
dependent surface tension parameterized in the follow- 
ing relation: cr(r) = cro[(T2 - T^)/{T^ + T2)]S/4,with 
(To = 18 MeV and T^ == 18 MeV (ct = at T > T^). The 
last contribution in Eq. Q involves the famous Fisher's 
term with dimensionless parameter r. 

The canonical partition function (CPF) of nuclear frag- 
ments in the SMM has the following form: 
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SiA-j:^knk). (2) 



In Eq. |2Il the nuclear fragments are treated as point- 
like objects. However, these fragments have non-zero 
proper volumes and they should not overlap in the co- 
ordinate space. In the excluded volume (Van der Waals) 
approximation this is achieved by substituting the to- 
tal volume V in Eq. (|2Jl by the free (available) volume 
Vf = V — h^i^kuk, where h = 1/pa [po = 0.16 fm~'^ 
is the normal nuclear density). Therefore, the corrected 
CPF becomes: Za{V,T) = Zf{V - bA,T). The SMM 
defined by Eq. |2Il was studied numerically in Refs. ,4;, 3|. 
This is a simplified version of the SMM, e.g. the symme- 
try and Coulomb contributions are neglected. However, 



its investigation appears to be of principal importance 
for studies of the liquid-gas phase transition. 

The calculation of Za{V,T) is difficult because of the 
constraint ^j, krik = A. This difficulty can be partly 
avoided by calculating the grand canonical partition 
function: 

oo 

Z{V,T,fi) = ^exp(i^)z^(t/,T)e(F-6A), (3) 

A=0 

where /i denotes a chemical potential. The calculation 
of Z is still rather difficult. The summation over {n^} 
sets in Za cannot be performed analytically because of 
additional A-dependence in the free volume V/ and the 
restriction Vj > 0. This problem was resolved [a, 13 
by the Laplace transformation method to the so-called 
isobaric ensemble lOj. 

In this work we would like to consider a more strict 

aV/b 

constraint ^ k nj- = A, where the size of the largest 

k 
fragment cannot exceed the total volume of the system 

(the parameter a is introduced for convenience). A sim- 
ilar restriction should be also applied to the upper limit 
of the product in all partitions Z]f(F,T), Zyi(y,T) and 
Z(V,T, /i) introduced above (how to deal with the real 
values of aV/b, see later). Then the model with this 
constraint, the CSMM, cannot be solved by the Laplace 
transform method, because the volume integrals can- 
not be evaluated due to a complicated functional In- 
dependence. However, the CSMM can be solved ana- 
lytically with the help of the following identity 
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which is based on the Fourier representation of the Dirac 
(5-function. The representation Q allows us to decouple 
the additional volume dependence and reduce it to the 
exponential one, which can be dealt by the usual Laplace 
transformation. Indeed, with the help of Q) the Laplace 
transform of the CSMM grand canonical partition (GCP) 
© can be done analytically: 



Z(A,T,^) = I dV e 
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After changing the integration variable V ^ V = V — 

a^/b 

b J2 k Uk, the constraint of 8-function has disappeared. 
k 



Then all Uk were summed independently leading to the 
exponential function. Now the integration over V' in 
Eq. |(SJ) can be straightforwardly done resulting in 
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where the function J-'{£,, A) is defined as follows 
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As usual, in order to find the GCP by the inverse 
Laplace transformation, it is necessary to study the struc- 
ture of singularities of the isobaric partition Q . 



III. ISOBARIC PARTITION SINGULARITIES 

The isobaric partition {Tj) of the CSMM is, of course, 
more complicated than its SMM analog [a, Q because for 
finite volumes the structure of singularities in the CSMM 
is much richer than in the SMM, and they match in the 
limit F ^ oo only. To see this let us first make the 
inverse Laplace transform: 
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dH^,x;j 



dXt. 



,(8) 



where the contour A-integral is reduced to the sum over 
the residues of all singular points A = A* 4- i?7 with 
n — 1,2, .., since this contour in complex A-plane obeys 
the inequality x > max(i?e{A*}). Now both remaining 
integrations in (|SJ| can be done, and the GCP becomes 






A* V 
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dT{V,Xl 



dXl 



(9) 



i.e. the double integral in JSJ simply reduces to the sub- 
stitution ^ — > V^ in the sum over singularities. This is a re- 
markable result which can be formulated as the following 
theorem: if the Laplace- Fourier image of the excluded vol- 
ume GCP exists, then for any additional V-dependence 
of !F{V,X'^) or 4>k{T) the GCP can be identically repre- 
sented by Eq. ^. 
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FIG. 1: A graphical solution of Eq. lO for T = 10 MeV and 
T = 1.825. The l.h.s. (straight line) and r.h.s. of Eq. (CU 
(all dashed curves) are shown as the function of dimensionless 
parameter I-i b for the three values of the largest fragment size 
aV/b. The intersection point at (0; 0) corresponds to a real 
root of Eq. llUt . Each tangent point with the straight line 
generates two complex roots of pOfl . 

The simple poles in ^ are defined by the equation 

a; = Hv,x*n). (10) 



In contrast to the usual SMM [6|, |3 the singularities A* 
are (i) functions of volume V, and (ii) they can have a 
non-zero imaginary part, but in this case there exist pairs 
of complex conjugate roots of H1Q(I because the GCP is 
real. 

Introducing the real i?„ and imaginary /„ parts of A* = 
Rn + ilii, we can rewrite Eq. (|10|l as a system of coupled 
transcendental equations (K^p ~ aV/b) 

Rn = J2 '^fe(^) ^^^^^ COs(/„6fc) , (11) 

fc=l 

/n = -^<^fe(r) e^^"^sin(/„6fc), (12) 

fc=i 

where we have introduced the effective chemical po- 
tential i^ = n + W{T) — \*ibT, and the reduced dis- 
tributions 4)i{T) = ziexp{-W{T)/T) and 0fe>i(T) = 
k~'^ exp(— cr(T) k'^'^/T) for convenience. 

Consider the real root (i?i > 0,/i = 0), first. For 
In ~ Ii — the real root i?i exists for any T and fi. Since 
for /„ 7^ defined by H12II the inequality cos(/„fefc) < 1 
cannot become the equality for all values of k simultane- 
ously, then from Eq. Hll|) one obtains 



Rn<J2^k{T) e- 



^ Rn < Rl 



(13) 



fc=i 



FIG. 2: Each curve separates the T — Reiy) region of one real 
root of Eq. IjlO^ (below the curve), three complex roots (on 
the curve) and five and more roots (above the curve) for three 
values of aV/b and the same parameters as in Fig. 1. 



(|13|l plays a decisive role in the thermodynamic limit 
y ^ oo because in this case it generates the pressure 
of the gaseous phase. 

Like in the usual SMM |a,l3, for infinite volume the ef- 
fective chemical potential must be real and non-positive, 
i/ < 0, because in the latter case the function J-{V, \\) ((Tj) 
diverges and the formal manipulations in |(SJ) to establish 
© cannot be used. The limiting value v = defines the 
liquid phase singularity of the isobaric partition which 
ives the liquid pressure pi{T, fi) = RiT == (^ -h W{T))/b 
Q|- But for finite volumes and finite Kup the effec- 
tive chemical potential can be complex (with either sign 
for its real part) and its value defines the number and 
position of the imaginary roots {A*>]^} in the complex 
plane. Consider the natural values of K^p, first. As it 
is seen from Fig. 1., the r.h.s. of (|12|) is the ampli- 
tude and frequency modulated sine-like function of di- 
mensionless parameter /„ b. Therefore, depending on T 
and Re{h') values there may exist either zero, or finite or 
infinite number of complex roots {A*>]^}. In Fig. 1. we 
showed a special case which corresponds to exactly three 



where the second inequality (|13|l immediately follows 
from the first one. Note that the second inequality 



roots of Eq. I|1U|) for each value of K^p'- the real root 
(/i = 0) and two complex conjugate roots (±/2)- Since 
the r.h.s. of (|12|) is monotonously increasing function of 
Re[u)^ when the former is positive, then it is possible to 
map the T — Reiy) plane into regions of a fixed number 
of roots of Eq. (|10|l . Each curve in Fig. 2. divides the 
T — Re{h') plane into three parts: for i?e(i')-values below 
the curve there is only a real root, for points on the curve 
there exist three roots, and above the curve there are five 
and more roots of Eq. (|10|) . 

A similar situation occurs for the real values of Kup- 
In this case all sums in Eqs. H1UI13|I should be expressed 



the Euler-MacLaurin formula 



fc=i { 



dkfik) 



/(A-„p)+/(2) 



A/(ifup)- 



A; (2), where A/(i^) = ^ 
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Here -B2n are the Bernoulh numbers. The representation 
(|14|l aUows one to study the effect of finite volume (FV) 
on the GCP ©. 



IV. FINITE VOLUME THERMODYNAMICS 

In the CSMM there are two different ways of how the 
finite volume affects thermodynamical functions: for fi- 
nite V and Reiy) there is always a finite number of simple 
poles in |j^ and all of them contribute into thermody- 
namic quantities; also the parameter a < 1 describes a 
difference between the geometrical shape of the volume 
under consideration and that one of the largest fragment 
(assumed to be spherical). To see this, let us study the 
mechanical pressure which corresponds to the GCP @ 



p = T 



dlTiZ{V,T,n) 
dV 



-Z(V,T,m) 
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J2MT)k'e'^+4>K^^{T) 
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'Kup [(1 - a) -f f (i - a)] -I- o{K^p) 



(15) 



where we give the main term for Aq and leading FV cor- 
rections explicitly for Re(y)/T < 1, whereas o{Kup) ac- 
cumulates the higher order corrections due to the Euler- 
MacLaurin Eq. (|14|l . In evaluation of H15|l we used an 
explicit representation of the derivative dX^/dV which 
can be found from Eqs. (|10|l and H14() . The first term in 
the r.h.s. of ((T3|) describes the partial pressure generated 
by the simple pole A* weighted with the "probability" 
e-^n "^ /Z{V, T, n), whereas the second and third terms ap- 
pear due to the volume dependence of K^p. Note that, 
instead of the FV corrections, the usage of natural values 



for Kupiy) would generate the artificial delta- function 
terms in (|15|l for the volume derivatives. 

As one can see from (|15f) for finite volumes the correc- 
tions can give a non-negligible contribution to the pres- 
sure because in this case Re{v) > can be positive. The 
real parts of the partial pressures TA* may have either 
sign. Therefore, according 113|) the positive pressures 
TRn > are metastable and the negative ones TRn < 
are mechanically unstable. The pair of complex conju- 
gate roots with the same value of TRn corresponds to a 
formation and decay of those states in thermodynamical 
system at finite volumes. 

When V increases the number of simple poles in © 
also increases and imaginary part of the closest to the 
real A-axis poles becomes very small. For infinite vol- 
ume the infinite number of simple poles moves toward 
the real A-axis to the vicinity of liquid phase singularity 
pi{T,ii)/T and, thus, generates an essential singularity 
of function T{V,pi/T) in jT)). In this case the contribu- 
tion of any of remote poles from the real A-axis to the 
GCP vanishes. Then it can be shown that the FV correc- 
tions in ((T^ become negligible because of the inequality 
Re[v) < 0, and, consequently, the reduced distribution 

of largest fragment 4'k^A'^) = ^up exp(-cr(r) Kli? /T) 
and the derivative dX^/dV vanish for all T-values, and 
we obtain the usual SMM solution [3, Q j and its thermo- 
dynamics is governed by the farthest right singularity in 
complex A-plane. The corrections of a similar kind should 
appear in the entropy, particle number and energy den- 
sity because of the T and /i dependence of A* due to 11U|) 
[ll|. Therefore, these corrections should be taking into 
account while analyzing the experimental yields of frag- 
ments. Then the phase diagram of the nuclear liquid-gas 
phase transition can be recovered from the experiments 
on finite systems (nuclei) with more confidence. 

Also it is possible that the metastable and unstable 
modes can emerge in the dynamically expanding system 
created in experiments, or even there may exist a direct 
relation between these modes and the spinodal instabil- 
ities discussed with respect to phase transition in finite 
systems [l^ , but it is theme of the other work. 
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